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ABSTRACT 

Aims. We report the discovery of CoRoT-8b, a dense small Saturn-class exoplanet that orbits a Kl dwarf in 6.2 days, and we derive its orbital 
parameters, mass, and radius. 

Methods. We analyzed two complementary data sets: the photometric transit curve of CoRoT-8b as measured by CoRoT and the radial velocity 
curve of CoRoT-8 as measured by the HARPS spectrometer**. 

Results. We find that CoRoT-8b is on a circular orbit with a semi-major axis of 0.063 ± 0.001 AU. It has a radius of 0.57 ± 0.02 Rj, a mass of 
0.22 ± 0.03 Mj, and therefore a mean density of 1.6 ± 0.1 gcirr 3 . 

Conclusions. With 67 % of the size of Saturn and 72 % of its mass, CoRoT-8b has a density comparable to that of Neptune (1.76 gcirr 3 ). We 
estimate its content in heavy elements to be 47-63 M e , and the mass of its hydrogen-helium envelope to be 7-23 M e . At 0.063 AU, the thermal 
loss of hydrogen of CoRoT-8b should be no more than ~ 0. 1 % over an assumed integrated lifetime of 3 Ga. 

Key words, planetary systems - stars: fundamental parameters - techniques: photometric - techniques: spectroscopic - techniques: radial veloc- 
ities 



1. Introduction 

CoRoT (Convection, Rotation and planetary Transits), a joint 
project of France, Austria, Belgium, Brazil, ESA, Germany, and 
Spain, is a low-Earth orbit visible photometer designed to mea- 
sure stellar light curves (LCs) for two main astrophysical pro- 
grams: asteroseismology and exoplanet detection (Baglin et al 



2006 ). Instrument characteristics and inflight performance can 
be found in Auvergne et al. (2009 1. CoRoT was launched on 
December 27, 2006 and will be operated by the French Centre 
National d'Etudes Spatiales (CNES) until March 31, 2013. 

The search for exoplanets in CoRoT LCs consists in look- 
ing for periodic dips caused by exoplanets whose edge-on orbits 
make them transit in front of their host stars. Study of transit 
signals combined with ground-based radial velocity follow-up 
observations leads to complete orbital solutions, as well as exo- 
planet masses and radii. It is with these systems that we learn the 
most about exoplanetary internal structures, atmospheres, and 
dynamics (for a review, see Charbonneau et al. 2007| l. Compared 



* Observations made with SOPHIE spectrograph at Observatoire de 
Haute Provence, France (PNP.07B.MOUT), and the HARPS spectro- 
graph at ESO La Silla Observatory (081.C-0388 and 083.C-0186). The 
CoRoT space mission, launched on December 27, 2006, has been de- 
veloped and is operated by the CNES with the contribution of Austria, 
Belgium, Brasil, ESA, Germany, and Spain. 

** Both data sets are available in electronic form at the CDS via 
anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
http : //cdsweb . u-strasbg . fr/cgi-bin/qcat? J/A+A/vol/page. 



to ground-based transit surveys, CoRoT offers uninterrupted 
photometric time series for up to about 150 days with relative 
errors as low as about 100 ppm per hour. This makes it possible 
to hunt for exoplanets that are only twice as big as the Earth in 
tight orbits, as demonstrated by the recent discovery of CoRoT- 
7b ( [Leger et al.|2009||Que"loz et al.|2009> . 

This paper reports the discovery and characterization of 
CoRoT-8b, a hot and dense small Saturn-class exoplanet in a 6.2- 
day orbit around a Kl dwarf star. We start by describing CoRoT 
LC analysis and transit parameter measurements. We then report 
ground-based follow-up observations and derive the full orbital 
and physical solution. After a look at transit timing variations 
and differential limbdarkening, we close this paper with a dis- 
cussion about CoRoT-8b's possible composition. 



2. CoRoT-8 photometric observations 

CoRoT performs photometric measurements by onboard sum- 
mation of all the photo-electrons collected during unit exposures 
of 32 s inside a predefined zone of the CCD matrix called a pho- 
tometric mask. At the beginning of a run, every target star re- 
ceives a mask whose shape is optimized given the star magnitude 
and environment. Mask sizes, ranging from about 4x4 to 15 x 10 
pixels or 9x9 to 23x35 arcsec on the sky, are fairly large in order 
to accomodate CoRoT defocused point-spread functions (PSFs). 
On the one hand, defocusing improves the photometric precision 
and reduces the sensitivity to telescope jitter, but on the other, it 
increases the probability of PSF overlaps that may bias the pho- 
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Table 1. Coordinates and magnitudes of CoRoT-8 and its six 
main contaminants. 



Star 


RA (deg) 


DE (deg) 


X 


Y 


V 


r' 


CoRoT-8 


291.58852 


1.42646 


924 


741 


14.8 


14.3 


Cont. #1 


291.58839 


1.42976 


929 


739 


18.8 


18.3 


Cont. #2 


291.58826 


1.42280 


919 


744 


19.1 


18.8 


Cont. #3 


291.58652 


1.42269 


920 


746 


20.7 


20.1 


Cont. #4 


291.59267 


1.42809 


924 


734 


20.7 


20.1 


Cont. #5 


291.58967 


1.43077 


929 


737 


19.0 


18.6 


Cont. #6 


291.58537 


1.42313 


921 


748 


20.3 


19.7 



Notes. X and Y are positions on CCD E2 in pixels. V and r' are the 
magnitudes in the Harris V and Sloan-Gunn r' filters. 



N ■ 



tometry (if left uncorrected) or cause transit false alarms (if a 
neighbor is an eclipsing binary). In the exoplanet channel, PSFs 
are not only defocused but also slightly dispersed by a biprism. 
Thus, for a subset of stars, CoRoT delivers three LCs instead of 
one by division of their masks into three submasks along the dis- 
persion direction. These "color LCs" are referred to as the red, 
green, and blue LCs, regardless of their true spectral content (and 
with no link whatsoever to any usual photometric system). 

CoRoT-8 is one of the 1 1 408 target stars observed by CoRoT 
during the first long run in constellation Aquila (LRcOl for 
short) from May 16, 2007 to October 5, 2007. |Cabrera et aT 
(2009} review all planetary transit candidates for this field, in- 
cluding CoRoT-8, which is referred to as LRcOl E2 1145-CHR- 
0101086161. This star belongs to the subset of 3719 stars for 
which color LCs are available (CHR stands for chromatic). The 
standard CoRoT pipeline (version 2.1 at this date) produces 4 
LCs in this case: one for every color and a "white" LC derived 
from the sum of colors. However, only the white LC receives 
proper telescope jitter noise correction, which is the reason color 
LCs have not been extensively studied in the previous papers of 
this series. For this study, we processed the white LC further 



by removing systematics effects as prescribed by Mazeh et al. 
(2009), and we applied a generalization of this treatment to the 
colors (Mazeh et al., in preparation). 

At this stage, none of these LCs is corrected for the fluxes of 
neighboring stars that leak into the target star mask. Exodat, the 
CoRoT entry database (Dele uil et al.| [2009), provides the nec- 
essary information to remove the contributions from these con- 
taminants. An estimate of the contamination for the white LC 
is readily available (1.6% in the case of CoRoT-8) but is com- 
puted with a generic mask prior to begining-of-run mask affecta- 
tion. To get a more precise estimate of the contamination for the 
white, as well as nonexisting estimates for the colors, we super- 
imposed inflight measured PSFs at the locations of the six main 
contaminants scaled with respect to their Sloan-Gunn r' mag- 
nitudes (Table [lj, then integrated their fluxes inside CoRoT-8's 
mask and submasks. Figure[T]displays the long exposure taken at 
the begining of LRcOl, CoRoT-8's mask, and the locations of the 
contaminants. For orientation, Fig. [2] is an image of the neigh- 
borhood of CoRoT-8 taken in 2004 with the Wide Field Camera 
(WFC) at the 2.5-m Isaac Newton Telescope (INT) at La Palma. 
Figure [3]is a computed version of Fig.[T|that ignores everything 
but CoRoT-8's and contaminants. This calculation yields new 
contamination estimates of 0.9 % for the white, 2.4 % for the 
blue, 0.2 % for the green, and 0.7 % for the red. In the following, 
LCs are corrected using these values. 

CoRoT-8's 142.1-day long LCs consist of a first section of 
11496 measurements with 512-s sampling period (sixteen 32-s 
unit exposures coadded onboard) and lasting 68.3 days. On the 



X-900 (pixel) 



Fig. 1. Long exposure taken at the beginning of LRcOl obser- 
vation showing CoRoT-8 and its environment. The mask cho- 
sen for CoRoT-8 (solid white line) contains a total of 55 pixels. 
The letters B, G, and R indicate the positions of the blue (20 
pix.), green (7 pix.), and red (28 pix.) submasks, respectively. 
The mask size on the sky is 1 8.6" x 16.2" (individual pixel size is 
2.32"). The white crosses indicate the positions of the six closest 
contaminants that lie inside a 40" -diameter circle (dashed white 
line) centered on CoRoT-8. X and Y are positions on CCD E2. 




-30 -20 -10 10 20 30 

RA- 291.58852 (arcsec) 



Fig. 2. Image taken in the r' filter with the WFC at the 2.5-m INT 
showing CoRoT-8 and its environment (seeing is ~ 2"). The six 
closest contaminants lying inside a 40" -diameter circle (dashed 
black line) are labeled 1 to 6 according to increasing angular 
distances from CoRoT-8. 



basis of this first section, the transit search pipeline called Alarm 
Mode ( Surace et al.|2008) l that continuously scrutinizes CoRoT 
LCs during a run, detected a transit-like signal and stopped 
unit-exposure onboard coaddition so that the remainder of this 
promising LC (198 752 measurements, 73.7 days) could bene- 
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X-900 (pixel) 



Fig. 3. Computed analog of Fig. [T] showing CoRoT-8's mask and 
its 6 closest contaminants using the same gray scale as Fig.[T] 



fit from a finer sampling period of 32 s. This action also trig- 
gered the beginning of ground-based follow-up operations: on- 
off transit photometry with the 1 .2-m telescope at Observatoire 
de Haute-Provence (OHP) in France showed that the transit sig- 
nal could not originate in any of the contaminants, so radial ve- 
locity (RV) spectrometry described in Sect. [5] was also started. 

Figure|4]shows the fully processed LCs (including systemat- 
ics removal and contamination correction) rebinned to a uniform 
512-s sampling period. The blue, green, and red fluxes amount 
to about 13 %, 14 %, and 73 %, respectively, of the total flux col- 
lected inside the photometric mask. Twenty-three transits with 
relative depths of about 1 % can easily be seen by eye in the 
white and the red, are more difficult to see in the green, and are 
barely visible in the blue. Discontinuities are caused by energetic 
particles (mostly protons) from the Van Allen radiation belts that 
hit CoRoT when it crosses the South-Atlantic Anomaly (SAA, 
Pinheiro da Silva et al.|20 08 ). They differ from one color to the 
next depending on the locations of the pixels impacted by SAA 
protons. The green LC seems to be more polluted than the other 
two. Overall, these LCs do not show obvious signs of stellar ac- 
tivity, in contrast to CoRoT-7, for instance, where quasi-periodic 
modulations could be used by Leger et al. ( 2009| l to estimate the 
stellar rotation period. 



3. White transit signal modeling 

We begin by measuring the transit signal properties by least- 
square fitting to the white LC of a periodic trapezoidal model 
with 5 parameters: the period P, the ingress start time T\, the 
ingress duration (taken as equal to the egress duration 734), 
the bottom-part duration 7*23, and the relative depth 8. 

We obtain first estimates of P, T\, the full transit duration 
T\4, and 5 by maximizing the cross-correlation with a periodic 
square-shape transit as explained by Borde et al. (2007). The 
white LC is then chopped in 23 transit-centered segments lasting 
3 X T\4. We normalize each segment and correct for local base- 
line flux variations by dividing by a parabola fitted to the out-of- 
transit flux variations. Then, we proceed to the least-square fit- 
ting of the 5-parameter trapezoid model described above, assum- 
ing measurement errors to be identical and equal to 8.7 10~ 4 , the 
standard deviation of out-of-transit measurements. We find that 
P = 6.2 1229 ±0.00003 days, T x = 2 454 238.9743+0.0004 HJD, 



6 = (6.51 ± 0.07) x lO" 3 , T n = 41.9 ± 0.8 min, and T 23 = 
80.4 ± 0.8 min with^ r 2 = 1.4. 

We keep the determined period and ingress start time for 
what follows, and we go on by fitting the white LC with a clas- 
sical model that describes the transit of an opaque black disk 
(the planet) in front of a limb-darkened (LD) luminous disk (the 
star). The normalized measured flux is / = 1 - F p /F+, where 
F p is the flux occulted by the planet and F* the flux emitted by 
the star. If I(x) is the normalized centro-symmetric LD profile of 
the star and x the normalized distance to the star center, then the 
total flux emitted by the star is = 2n J Q I(x) x dx. Denoting 
r = Rp/R* as the ratio of the planet radius to the star radius, the 
flux F p occulted by the planet at linear distance p from the star 
center can be computed via 



r p+r (x 2 +p 2 -r 2 \ 

2 I I(x) arccos x dx 

Jp-r \ 2px / 

if r < p < 1 + r, and with 

J^r-p rr+p I , 

I(x) xdx + 2 I I(x) arccos - 
Jr-p \ 



(1) 



' x 2 + p 2 — r 2 \ 

2n I I(x) xdx + 2 I I(x) arccos I x dx (2) 



2px 



if < p < r. Besides, 



7 ( 2n(t-T Q ) \ | . 2 ( 2n(t-T ) \ 
9 = R~ V 1 p J ( p )' ( 

where a is the semi-major axis, i the orbit inclination, P the or- 
bital period, and Tq the central time of the first transit. To con- 
sider the effect of signal integration over 512-s exposures, our 
model is computed every 8 s, then numerically integrated over 
64 points. In the course of our work on CoRoT planets, this 
set of equations and its numerical implementation have been 
checked against other published approaches (Mand eT& Agol| 
[20021 |GimenezJ2006| ). 

From Eqs.[TJ|3j we see that this model's parameters are the 
radius ratio R p /Ri,, the scaled semi-major axis a/R+, the orbital 
inclination i, the transit central time r , and the orbital period P, 
plus usually \-A LD parameters depending on the form given to 
I(x). In the following, we do not fit Tq and P, but instead take the 
values obtained with the trapezoid model, and use the quadratic 
LD law I(p) = 1 - mi (1 - yu) - u 2 (1 - /j) 2 (where the fitted LD 
parameters are actually u + = u\ +U2 and m_ = u\ Therefore, 
we have a total of 5 parameters to determine. 

Parameter optimization, done by way of a genetic algorithm 

leads to R P IR+ = 

0.075 ± 0.001 



a/R* 



88.4 



±0.1°, 



7.6 ± 0.4, ; 

0.09, and u 2 = 0.14 + 0.09 with*, 2 = 1.4 (Fig. [5} 
parameter errors were found by increasing^- 2 from^ 



u x = 0.70 ± 
Individual 



min 



toy 



1 as prescribed by [Press et aL] ( |1997| l. From the scaled semi- 
major axis and the period, we obtain (via Kepler's third law) 



1 /3 

MJ /R* = 1.24 ± 0.03 in solar units, hence a stellar density 
p-i, = 2.7 + 0.1 g cirT 3 (see Table[5]for a summary). To compute 
the planet radius, we now need to discuss the star properties. 



4. C0R0T-8 stellar parameters 

C0R0T-8 was observed with the UVES spectrograph on the VLT 
(program 081.C-0413) on September 2 and 3, 2008. We used 
the Dicl mode (346+580) and a slit width of 0.5" that pro- 
vides a resolving power of ~ 65 000. The total exposure time 
was 4760 s, divided into two exposures. The data were reduced 
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Fig. 4. CoRoT-8's normalized light curves vs. time in Heliocentric Julian Date (HJD) minus 2451 545 (January 1 st , 2000). Sampling 
period is 512 s. The actual median fluxes are 14724, 15 677, 80528, and 110 408 e~/32 s for the blue, green, red, and white LCs, 
respectively. 
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Fig. 5. CoRoT-8b's best transit model for the white LC (solid 
line) with individual measurements (small dots) and averaged 
measurements inside 10-min bins (big dots). Fit residuals are 
displayed around 0.985. 



using our own reduction software to optimize the spectrum ex- 
traction. The spectra were divided by the blaze function and put 
at rest by subtracting the systemic radial velocity. We obtained a 



signal-to-noise ratio (S IN) of about 170 per resolution element 
at 550 nm. 

The fundamental parameters of the star were first estimated 
with this spectrum. Later, we took advantage of the series of 
HARPS spectra to refine the estimates of v sin i and the micro 
and macro-turbulent velocities, which all benefit from the much 
higher spectral resolution of HARPS. To that purpose, among 
the HARPS spectra collected for the radial velocity analysis 
(Sect. |5]l, we selected those that were not contaminated by Moon 
light. Nine spectra were finally co-added, yielding an average 
spectrum with S /N ~ 150 per resolution element at 550 nm. 

The star is a slow rotator since the rotational broadening is 
lower than the spectral resolution of HARPS. Using the synthe- 
sis of a set of isolated lines, we determined a projected rotational 
velocity of v sin i = 2 ± 1 km s , and a macroturbulent velocity 
of v mac = 2+1 kms _I . For the spectroscopic analysis, we fol- 
lowed the same methodology as for the previous CoRoT planets 
since CoRoT-3b. The analysis was performed using two spectral 
synthesis packages VWA (|Bruntt et al.|2002) and SME ( |Valenti 
|& Piskunov||1996 |Valenti & Fischerp 005) as described in de- 
tails by |Bruntt| ( |2009| l. 

For such a late type star that is metal rich, the continuum 
placement could be a serious issue at wavelengths shorter than 
550 nm owing to molecular bands and the high level of blend- 
ing of weak metal lines. We therefore concentrated our spectro- 
scopic analysis on well-defined spectral lines at longer wave- 
lengths. We used the Cai lines at 612.2, 616.2, and 643.9 nm 
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Table 2. Abundances of the main elements in CoRoT-8. 



CoRoT-Bb 



HARPS + SOPHIE 



Element 


[X/H] 


number 






of lines 


Mgi 


0.16 (0.16) 


3 


Ca i 


All /A AC\ 

0.31 (0.05) 


4 


Til 


0.26 (0.05) 


14 


V i 


0.49 (0.07) 


16 


Fei 


0.31 (0.05) 


96 


Fe ii 


0.32 (0.06) 


8 


Co i 


0.30 (0.07) 


6 


Nil 


0.34 (0.05) 


30 


Si i 


0.37 (0.05) 


18 



as surface gravity diagnostics instead of the usual Mgib or Nai 
lines. This overcomes the issue of continuum placement for such 
broad lines with extended wings when observed with an echelle 
spectrum. To provide an additional check of our effective tem- 
perature determination, we compared the Ha line profile outside 
of the core to theoretical profiles, following the method recom- 
mended in Sect. 3.3 of Bruntt et al. ( 2004| l. The result agrees with 
our estimates obtained independently with VWA and SME. Our 
final atmospheric parameters, T e ff , log g, and global metallicity 
estimates are given in Table [5] The results from the detailed el- 
emental abundances analysis are summarized in Table [2] These 
correspond to a slowly rotating Kl main-sequence star with no- 
ticeable metal enrichment. 

To estimate the mass and the radius of the host star, we com- 
bined the constraints from the spectroscopic parameters (r e jf and 
[M/H]) and M l J 3 /R+, the proxy for stellar density derived from 
the transit physical model (see Sect. [3] and Table [5J. The spec- 
troscopic analysis of our spectra gave no evidence of youth: the 
Li i line at 670.78 nm is absent and the star's spectra display no 
hint of chromospheric activity. This is consistent with the low ro- 
tational velocity we measured. Using STAREVOL stellar evolu- 
tionary tracks ( Siess 2006 , Palacios, private communication), we 
thus excluded evolutionary tracks corresponding to pre-main se- 
quence stages and assumed that the star is on the main sequence. 
In this way, we obtained a stellar mass of M* = 0.88 ± 0.04 M Q 
and furthermore inferred a stellar radius of R+ = 0.77 + 0.02 R . 
The corresponding surface gravity is 4.61 + 0.07 in good agree- 
ment with the spectroscopic value of log g = 4.58 + 0.08. 

Concerning limb-darkening, we also note fair agreement be- 
tween measured values reported in Sect. [3] u\ = 0.70 ± 0.09 
and t<2 = 0.14 ± 0.09, and theoretical estimates computed for the 
CoRoT bandpasses by Sing (2010): u\ = 0.59 and u\ = 0.12 (for 
r eff = 5000 K, logg = 4.5, and [M/H]=0.3). Achieving consis- 
tency in both surface gravity and limb-darkening is a good in- 
dication that the determined orbital and physical parameters are 
reliable. Therefore, we can safely derive a planetary radius of 
0.57 + 0.02 Rj, which is 1.7 times the radius of Neptune and 0.7 
times that of Saturn. 



5. Radial velocity spectrometry 

We performed RV observations of CoRoT-8 with the SOPHIE 
spectrograph (Perruchot et al. 2008 1 Bouchy et al. 2009a| l at 
the 1.93-m telescope located at OHP, and the HARPS spec- 
trograph (Mayor et al. |2003 1 at the 3.6-m ESO telescope at 
La Silla Observatory (Chile). These two instruments are cross- 
dispersed, fiber-fed, echelle spectrographs dedicated to high- 
precision Fizeau-Doppler measurements. SOPHIE was used in 
its high-efficiency mode (spectral resolution of 40000). HARPS 
and SOPHIE were both used in the obj_AB observing mode, 



-37.4 



,— , -37,45 
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Fig. 6. Period-folded radial velocities of CoRoT-8. The black 
dots and open red squares respectively correspond to HARPS 
and SOPHIE measurements. 



where the second fiber was made available for monitoring the 
background Moon light at the cost of no simultaneous Th-Ar 
lamp calibration. This is made possible thanks to the remark- 
able intrinsic stability of these spectrographs. Indeed, the con- 
tribution of instrumental drift to the global RV errors was al- 
ways negligible with respect to the photon noise error. HARPS 
and SOPHIE data were reduced with the standard pipeline based 
on the work by Baranne et al. ( 1996 1, Pepe et al. (2002), and 
Bouchy et al. (2009a I. In the case of CoRoT-8, radial veloci- 
ties were derived by cross-correlating the acquired spectra with 
a numerical K5 mask. 

The first two measurements on CoRoT-8 were made with 
SOPHIE on August 2007 and showed no significant RV vari- 
ations at the level of 50 ms _1 . As such, they were compatible 
with a planetary companion lighter than 1 Mj, so it called for 
more observations. Nine more measurements were obtained with 
HARPS between the end of July 2008 and the end of September 
2008 (ESO program 081.C-0388). Although they yielded RV 
variations in phase with the ephemerides computed from the 
CoRoT-8 light curve, they did not allow by themselves properly 
constraining of the orbit, and above all completely ruling out the 
presence of a blended binary (a background eclipsing binary not 
angularly separated from CoRoT-8 by the 3.6-m ESO telescope). 
Therefore, ten additional HARPS measurements were acquired 
to secure the planetary origin of the CoRoT transits from June to 
August 2009 (ESO program 083.C-0186). 

All RV measurements (Table|3]l are displayed in Fig. [6] The 
SOPHIE data were shifted by 105 ms 'for this plot. Several 
spectra of the third campaign were slightly contaminated by 
Moon light. We attempted to correct for the Moon signal, but 
it only introduced more noise and did not significantly affect the 
target RV values. (The Moon RV was different enough from that 
of the target.) We concluded from this study that Moon contam- 
ination could be safely ignored in our case. 

The period-folded radial velocities in Fig. [6] present clear 
variations in phase with ephemerides derived from CoRoT tran- 
sits, in agreement with what would be expected for the reflex 
motion due to a planetary companion. We fitted RV data with 
a Keplerian model assuming a circular orbit and found a sys- 
temic RV of -37.468 ± 0.003 km s 1 and an RV semi-amplitude 
of 26 ± 4 ms^'with^ = 0.8 (showing that the standard devi- 
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ation of the residuals equal to 11.8 ms 'has the same order of 
magnitude as the individual RV errors). 

To examine the possibility that the radial velocity variation 
is caused by a blended binary, we followed the procedure de- 
scribed in Bouchy et al. (2009b). It consists in checking the 
spectral line assymmetries, as well as the dependence of RV 
variations on the cross-correlation mask used to compute the 
cross-correlation function (CCF). Spectral line assymetries are 
revealed by the calculation of the bisector span of the CCF. To 
increase the signal-to-noise ratio, we first averaged CCFs around 
orbital phases of 0.25, 0.5, and 0.75: six, five, and five CCFs 
were thus averaged in orbital phase intervals 0.15-0.35, 0.40- 
0.61, and 0.65-0.86, respectively. Bisector spans computed on 
these three averaged CCFs show no significant correlation with 
the RV as can be seen in Fig. [7] since the rms value of the bi- 
sector span variation is less by a factor of five than the radial- 
velocity semi-amplitude K. Furthermore, RV computed with dif- 
ferent cross-correlation templates (for F0, G2, and K5 stars) do 
not correlate with the amplitude of RV. 

These two checks are common tools for discarding the 
blended binary scenario and confirming the planetary nature of 
transiting planets. However, to support this conclusion further, 
we also contemplated the two scenarios: 1) a triple system, and 
2) a background eclipsing binary blended with CoRoT-8. In the 
first case, it was found that only a brown-dwarf companion or- 
biting an M dwarf could produce the same RV and bisector sig- 
nals. However, such a system would have a relative transit depth 
twice as large as the one measured by CoRoT. In the second case, 
the velocity of CoRoT-8 and the systemic velocity of the back- 
ground binary would have to be the same at the level of 1 km/s 
in order to give a symmetrical spectroscopic signal. This would 
be very unlikely when considering the dispersion of the velocity 
distribution in the Galaxy. Therefore, the most probable cause of 
the RV signal is a planetary companion to CoRoT-8. 

CoRoT-8b's mass (in Mj) can then be computed as 
4.92 10~ 3 KM 2 J 3 P l 13 , where K is in ms" 1 , M* in M , and P 
in days. We find M p = 0.22 + 0.03 Mj, 4.1 times the mass of 
Neptune and 0.7 times that of Saturn. Other stellar and planetary 
parameters can be found in Table [5] 

Finally, no significant offset was found between the two sets 
of HARPS measurements separated by one year, thus excluding 
the presence of an additional massive giant planet at a few AUs 
in the system. 



6. Search for another body by transit timing 

In this section, we look for transit timing variations (TTVs) in 
the white LC that may be caused by a gravitational perturber 
(e.g., |Sartoretti & Sclmeider|1999||Schneider|2003| >. To do this, 
we fit a trapezoid to each transit separately in which the ingress 
time, T\ , is the only free parameter and other parameters are held 
constant at the values reported at the beginning of Sect [3] 

We find variations that seem to be statistically signifi- 
cant (Fig. [8} and whose autocorrelation curve peaks at IP = 
43.5 days. However, we cannot be assured that these variations 
are gravitational in origin. To see whether they might be caused 
by stellar activity in the LC, we performed a dedicated analy- 
sis of the LC long-term variations to measure the stellar rota- 
tion period. For this purpose, we selected the red LC as the best 
compromise between the photometric precision and the number 
of discontinuities (Fig. [4]). We corrected this LC for two main 
discontinuities around dates 2733 and 2824 and removed data 
points inside the transits. The autocorrelation of the resulting LC 
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Fig. 7. Bisector span vs. radial velocity. The three points shown 
are averages obtained from six, five, and five cross-correlation 
functions computed in orbital phase intervals 0.15-0.35, 0.40- 
0.61, and 0.65-0.86, respectively. 

Table 3. Radial velocity measurements on CoRoT-8 obtained by 
HARPS and SOPHIE spectrometers. 



BJD° 


RV 


1-cr error 


exp. time 


S/N p. pix. 


-2400000 


(km s" 1 ) 


(kms- 1 ) 


(s) 


(at 550 nm) 


SOPHIE 


54314.40776 


-37.6060 


0.0490 


3600 


12.15 


54318.41015 


-37.5454 


0.0184 


3600 


17.00 


HARPS 


54675.69652 


-37.5050 


0.0151 


2700 


9.20 


54678.66557 


-37.4427 


0.0101 


2700 


12.20 


54679.73180 


-37.4751 


0.0407 


2700 


4.60 


54682.68436 


-37.4830 


0.0112 


2700 


11.20 


54700.57459 


-37.4756 


0.0140 


2700 


8.90 


54703.59565 


-37.4611 


0.0169 


2700 


8.40 


54732.56202 


-37.4687 


0.0080 


3600 


14.50 


54733.58576 


-37.4508 


0.0077 


3600 


14.80 


54737.58497 


-37.4955 


0.0134 


3600 


10.30 


54987.78368 


-37.4550 


0.0237 


3600 


7.82 


54988.84456 


-37.4482 


0.0146 


3600 


10.92 


54989.83528 


-37.4540 


0.0098 


3600 


14.77 


54991.85115 


-37.4979 


0.0157 


3600 


10.74 


54992.80807 


-37.4839 


0.0106 


3600 


13.68 


54993.87346 


-37.4610 


0.0113 


3600 


13.30 


54998.85102 


-37.5006 


0.0106 


3600 


12.20 


55021.82951 


-37.4656 


0.0102 


3600 


14.22 


55045.74089 


-37.4139 


0.0148 


3600 


10.95 


55046.73467 


-37.4888 


0.0179 


3600 


9.47 


(a) Barycentric Julian Date 



has a weak local maximum around 20 days (20 ± 5 days). If we 
cut this 142-day LC into three thirds lasting about 47 days each, 
we again find weak local maxima around 20 days for the first and 
the third thirds, but not for the second one. Although uncertain, 
this tentative signal around 20 days might still reflect the rota- 
tion period of the star. Indeed, if we combine the measurements 
for v sin R+, and i, we derive a loosely constrained stellar rota- 
tion period of 20 ± 10 days. Now, activity-induced TTVs may be 
periodic if the planet transits repeatedly above the same stellar 
surface features, which implies that the ratio of the stellar rota- 
tion period to the planet orbital period be a rational number. If 
7P/2 = 21.7 days is closer to the true stellar rotation period than 



6 



P. Borde et aL: Transiting exoplanets from the CoRoT space mission 



Transit timing variations 



-1 




Table 4. Comparison of limb-darkening parameters measured 
from the white and color LCs. 
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Fig. 8. Transit-to-transit variation of the date of first contact as a 
function of the transit number. 
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Fig. 9. Transit-to-transit variations in the relative depth, date of 
first contact, and first to fourth contact duration as a function of 
the transit number. 



3P = 18.6 days, this may explain why there seem to be a TTV 
signal at about twice the stellar rotation period. 

We also explored transit-to-transit variations of 6, T\, and 
T\4 at the same time by fitting trapezoids to individual transits 
where 8, T\, T\i, and T23 are now all free parameters. Figure [9] 
displays the resulting variations, denoted A6, AT\ , and ATu, as 
a function of the transit number. The last two variations are ob- 
viously highly anti-correlated (with a correlation coefficient of 
-0.82) showing that variations in date actually reflect variations 
in duration. As this should not be the case for a gravitational per- 
tuber, we are inclined to believe that the TTVs we measure are 
more likely to be the effect of stellar noise. 



7. Differential limb-darkening 

With its photometry in three colors, CoRoT makes it possible in 
principle to study differential stellar limb-darkening. However, if 
the photometric accuracy in the red is comparable to the white, 
it is not the case for the green and the blue where the col- 
lected flux is about five times lower than in the red (Sect. [2j. 



LC 


U\ 


«2 




white 


0.70 ± 0.09 


0.14± 


0.09 


blue 


0.88 ± 0.04 


-0.04 i 


:0.04 


green 


0.68 ±0.16 


-0.25 i 


:0.16 


red 


0.46 ±0.12 


0.26 ± 


0.12 



Indeed, the standard deviation of out-of-transit measurements 
is 8.0 10~ 4 for the red, and jumps to 2.4 10~ 3 for the blue and 
2.1 10" 3 for the green. Moreover, satellite jitter noise and con- 
tamination by neighboring stars (Sect. [2J increase the difficulty 
of this endeavor. Nevertheless, we make a first attempt here to 
study CoRoT-8b's transit in the CoRoT colors. 

As a start, we first used the same trapezoidal model as for the 
white but we fixed P, T\, and Tu at the values obtained with the 
white LC (as we do not expect variations here), and fit only for 6 
and T\i- We find that the transit ingress duration is longer in the 
blue (5 1 min) than in the green (46 min) and longer in the green 
than in the red (39 min), which is the expected effect of the more 
pronounced star limb-darkening at shorter wavelengths, and this 
seems encouraging. The blue transit ((5.71 +0.17) x 10~ 3 ) is also 
slightly deeper than the green ((5. 14 + 0.14) x 10~ 3 ), as expected, 
but is unfortunately not deeper than the red ((5.71 ±0.17)x 10~ 3 ), 
which suggests that color corrections might still not be accurate 
enough for CoRoT-8. 

Nevertheless, if we attempt to model the color transits with 
the white model from Sect.[3]by holding constant R p /R+, a/R+, 
and i, and letting u\ and «2 vary, we do measure limb-darkening 
parameters for the colors (Table 0), but with a moderate fit qual- 
ity in the green and the blue. 

It also worth reporting that, when R p /R+, a/R*, and i are 
included as free parameters in the modeling of the three colors 
taken separately, the obtained best-fit values for these parame- 
ters are all consistent with the white solution. This consistency 
check brings further evidence that the transiting object is a gen- 
uine exoplanet. 



8. Discussion 

With a transit depth of about 7 mmag, a transiting planet like 
CoRoT-8b is almost impossible to detect for ground-based pho- 
tometric surveys that are not sensitive to giant planets smaller 
than Saturn: only HAT-P-1 1 b ( |Bakos et al.|2010j ) with a smaller 
transit depth (4.3 mmag) was detected by a ground-based pho- 



tometric survey. CoR oT-7b (0.3 mmag, |Leger et al.||2009| l and 
Kepler-4b (0.9 mmag, Boruc ki et al.|2010| l were both detected by 
space-based photometric surveys, whereas GJ 436 b (0.6 mmag, 
Gillon etal.|2007) and HD 149026 b (2.6 mmag, |Satoetal.|2005) 
were both detected by RV surveys prior to their photometric ob- 
servations. 

With the exception of CoRoT-10 (V = 15.2, Bonomo et al., 
in preparation), CoRoT-8 (V = 14.8) is the faintest star harboring 
a planet found by the CoRoT survey. It is also relatively faint by 
RV standards, so that about 20 HARPS measurements (totaliz- 
ing about 20 hours of telescope time) split over two campaigns 
were needed to definitively establish the planetary nature and 
constrain the mass of CoRoT-8b within 14 %. 

Radial velocity surveys do not find a lot of planetary com- 
panions at close distances (P < 10 days) with masses compara- 
ble to CoRoT-8b. In the mass-period diagram (Fig. 10 1, CoRoT- 
8b appears in between hot-Jupiter and super-Earth populations. 
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Table 5. Star and planet characteristics of the CoRoT-8 system. 



CoRoT-ID 
CoRoT-WinID 

GSC 2.3 

Coordinates (J2000) 

Magnitudes B, V, r', i' 

Results from light-curve analysis 

Planet period, P 
Transit epoch, To 
Transit duration, T l4 
Transit relative depth, 6 
Radius ratio, Rp/R* 
Scaled semi-major axis, a/R* 
Orbital inclination, i 
Impact parameter, b 
Limb-darkening coefficients, u\, ui 
Stellar density, p+ 
M'f/Rt 



0101086161 
LRc01-E2-1145 

N29S083644 

19:26:21.245 +01:25:35.55 
16.10, 14.80, 14.27, 13.41 

6.21229 ± 0.00003 days 
HJD 2 454 238.9743 ± 0.0004 
2.74 ± 0.02 h 
(6.51 ± 0.07) x 10- 3 
0.075 ± 0.001 
17.6 ±0.4 
88.4 ±0.1° 
0.49 ± 0.04 

0.70 ±0.09, 0.14 ±0.09 
2.7 ±0.1 gem- 3 



1.24 ±0.03 (solar units) 
Results from radial velocity observations 



26 ± 4 ms-' 
(fixed) 

-37.468 ±0.003 km s-' 



Radial velocity semi-amplitude, K 
Orbital eccentricity, e 
Systemic radial velocity 
Results from spectral typing 

Spectral type 

Effective temperature, T^g 
Surface gravity, log g 
Metallicity, [M/H] 
Micro-turbulent velocity, v mic 
Macro-turbulent velocity, v mac 
Rotational velocity, v sin i 
Star mass, M t 
Star radius, R+ 
Star age 
Star distance 

Absolute physical parameters from combined analysis 

Planet mass, M p 0.22 ± 0.03 M, 

Planet radius, R p 0.57 ± 0.02 Rj 

Planet density, p p 1 .6 ± 0. 1 g cm 3 

Planet orbital semi-major axis, a 0.063 ± 0.001 AU 



This could challenge the existence of a bi-modal distribution 
(Mordasini et al. 2009), or CoRoT-8b might belong to the dis- 
tribution tail of giant gaseous planets. 

CoRoT-8b has a density of 1.6 gem- 
higher than for Saturn (0.69 g cirT 3 ), but comparable to Neptune 
(1.76 gem 3 ) or HD 149026b (1.7 gcirT 3 ). The last orbits at 



K1V 

5080 ± 80 K 
4.58 ± 0.08 (cgs) 
0.3 ±0.1 
0.8 ± 0.2 km/s 
2 ± 1 km/s 
2± 1 kms 1 
0.88 ± 0.04 M G 
0.77 ± 0.02 R Q 
<3Ga 
380 ± 30 pc 



(Fig. 



11 



which is 



0.043 AU from its host star and has similar mass and radius to 
CoRoT-8b (0.36 Mj, 0.65 Rj). By contrast, the density of HAT- 
P-12b (0.21 Mj, 0.96 Rj) that orbits at 0.038 AU from its host 
star is only 0.32 gcirT 3 . 

We estimate the thermal loss of hydrogen from CoRoT-8b 
by applying the method described by Lammer et al. (2009). 
According to the expected evolution of the soft X-rays and EUV 
flux of K stars and CoRoT-8b's orbit location at 0.063 AU, the 
planet lost not more than about 0.045-0.12 % over an assumed 
integrated lifetime of 3 Ga. The lower mass-loss value corre- 
sponds to a heating efficiency of 10 % and the higher value to a 
heating efficiency of 25 %. From this estimation we can conclude 
that CoRoT-8b has only slightly evolved since its formation. 

To constrain the planetary composition, we combine stellar 
and planetary evolution models. The stellar evolution models are 
based on CESAM (Morel & Lebreton|2008] l, include the diffu- 
sion of chemical elements in the star's radiative zone, and are 
calibrated to match our Sun for a mass of 1M , a metallicity 
[M/H] = 0, and an age of 4.5 Ga. A grid of models is com- 
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Fig. 10. Correlation diagram between mass and period for close- 
in planets. CoRoT-8b appears to lie in between the bulk of giant 
planets and the bulk of super-Earths. Open triangles correspond 
to planets with unknown orbital inclinations. 
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Fig. 11. Mass-radius diagram for transiting planets including 
CoRoT-8b. The two dashed lines correspond to the densities of 
Saturn (0.69 gcrrT 3 ) and Jupiter (1.33 gem- 3 ). Neptune has a 
density of 1.76 gcirT 3 . 

puted as a function of mass, metallicity, and age. The two major 
constraints used are CoRoT-8's stellar density and effective tem- 
perature. For a given model, we define 



r eff (M*, [M/H], (J y + /p*(M*, [M/H], O 

°~t cB 



1/2 



(4) 



which may be considered as the number of standard devia- 
tions from the combined (r e ff,p*) constraints with correspond- 
ing standard deviations (o-T cfr ,<T Pt ). We look for models where 
yia- t is small and at least less than 3. (The constraint on the 
spectroscopically-determined logg value was found to be too 
weak to be useful.) Among these models, we select only the 
ones which have the right metallicity at the age considered, i.e., 
0.2 < [M/X](f) < 0.4. Since we account for chemical diffusion, 
the value of [M/X](f) gets progressively lower with time. For a 
metal-rich star like CoRoT-8b, this favors slightly younger ages, 
but the star's low mass makes the effect very small. 

The corresponding constraints on the planetary parameters 
are then derived using the photometric and radial velocity con- 
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Fig. 12. Evolution of the size of CoRoT-8b (in Jupiter units, 
1 Rj= 71492 km) as a function of age (in Ga), compared 
to constraints inferred from CoRoT photometry, spectroscopy, 
radial velocimetry, and stellar evolution models. Red, blue, 
and green dots correspond to the planetary radii and ages 
that result from stellar evolution models matching the inferred 
(p+, r e ff)-uncertainty ellipse within ltr, 2<x, and 3cr, respectively. 
Planetary evolution models for a planet with a solar-composition 
envelope over a central dense core of pure rocks and pure ices of 
variable mass are shown as plain and dashed lines, respectively. 
These models assume an equilibrium temperature of 925 K and 
a total mass of 0.22 Mj. The inferred solutions at a lcr level are 
bounded by models with a minimum (rock) core mass of 47 M ffi 
and a maximum (ice) core mass of 62.5 M®. The corresponding 
masses of the hydrogen-helium envelope are 22.9 M ffi (32.8 % 
M tot ), and 7.4 M® (10.6 % M tot ), respectively. 



straints. Planetary evolution models are calculated with standard 
hypotheses (e.g., |Guillot|[2008[ ): we assume that the planet is 
made of a central rock/ice core of variable mass and of an over- 
laying envelope of solar composition. Given that the planet is 
farther from its star than most other transiting planets, we do not 
consider the possibility of additional energy input from stellar 
tides. 

The results in terms of planetary size as a function of sys- 
tem age are shown in Fig. [12] The colored areas indicate the 
constraints derived from the stellar evolution models, for values 
of «o- 4 = 1,2, and 3, respectively. The age constraints are rel- 
atively weak but favor values below 3 Ga. The planetary radius 
is tightly constrained, on the other hand. The amount of heavy 
elements that are found to bound the = 1 surface are 47 to 
56 M ffi of rocks or 56 to 62.5 M ffi of ices. 

With a density comparable to that of Neptune, CoRoT-8b has 
a mass of heavy elements similar to that of HD 149026 b (Sato 
eTaT1[2005| |Ikoma et al.||2006 



but a much smaller hydrogen- 
helium envelope. As for HD 149026, the parent star is metal- 
rich, confirming the trend that metal-rich stars tend to form gi- 



ant planets with high contents in heavy elements ( |Guillot et al. 
2006, |Burrows et aT]|2007l |Guillot||250"8] l. The question of the 
formation of these planets with high masses in heavy elements 
but small hydrogen-helium envelopes remains open, but proba- 
bly involves giant collisions ( Ikomaet al.|2006 i. 



Acknowledgements. We wish to thank the French National Research Agency 
(ANR-08-JCJC-0 102-01) for its continuous support for our planet-search pro- 
gram. The team at the IAC acknowledges support by grant ESP2007-65480-C02- 
02 of the Spanish Ministerio de Ciencia e Innovacion. The German CoRoT Team 



References 

Auvergne, M., Bodin, P., Boisnard, L., et al. 2009, A&A, 506, 41 1 

Baglin, A., Auvergne, M., Barge, P., et al. 2006, in ESA Special Publication, 

Vol. 1306, ESA Special Publication, ed. M. Fridlund, A. Baglin, J. Lochard, 

& L. Conroy, 33-+ 
Bakos, G. A., Torres, G, Pal, A., et al. 2010, ApJ, 710, 1724 
Baranne, A., Queloz, D., Mayor, M., et al. 1996, A&AS, 1 19, 373 
Borde, P., Fressin, F, Ollivier, M., Leger, A., & Rouan, D. 2007, in Astronomical 

Society of the Pacific Conference Series, Vol. 366, Transiting Extrapolar 

Planets Workshop, ed. C. Afonso, D. Weldrake, & T. Henning, 145 — h 
Borucki, W. J., Koch, D. G, Brown, T. M., et al. 2010, ArXiv e-prints 
Bouchy, E, Hebrard, G, Udry, S., et al. 2009a, A&A, 505, 853 
Bouchy, E, Moutou, C, Queloz, D., & the CoRoT Exoplanet Science Team. 

2009b, in IAU Symposium, Vol. 253, IAU Symposium, 129-139 
Bruntt, H. 2009, A&A, 506, 235 

Bruntt, H., Bikmaev, I. E, Catala, C, et al. 2004, A&A, 425, 683 

Bruntt, H., Catala, C, Garrido, R., et al. 2002, A&A, 389, 345 

Burrows, A., Hubeny, I., Budaj, J., & Hubbard, W. B. 2007, ApJ, 661, 502 

Cabrera, J., Fridlund, M., Ollivier, M., et al. 2009, A&A, 506, 501 

Charbonneau, D., Brown, T. M., Burrows, A., & Laughlin, G. 2007, Protostars 

and Planets V, 701 
Deleuil, M., Meunier, J. C, Moutou, C, et al. 2009, AJ, 138, 649 
Geem, Z., Kim, J., & Lonatathan, G. 2001, Simulation, 76, 60 
Gillon, M., Pont, E, Demory, B., et al. 2007, A&A, 472, L13 
Gimenez, A. 2006, A&A, 450, 1231 
Guillot, T. 2008, Physica Scripta Volume T, 130, 014023 
Guillot, T., Santos, N. C, Pont, E, et al. 2006, A&A, 453, L21 
Ikoma, M., Guillot, T., Genda, H., Tanigawa, T., & Ida, S. 2006, ApJ, 650, 1 150 
Lammer, H., Odert, P., Leitzinger, M., et al. 2009, A&A, 506, 399 
Leger, A., Rouan, D., Schneider, J., et al. 2009, A&A, 506, 287 
Mandel, K. & Agol, E. 2002, ApJ, 580, L171 

Mayor, M., Pepe, E, Queloz, D., et al. 2003, The Messenger, 114, 20 

Mazeh, T., Guterman, P., Aigrain, S., et al. 2009, A&A, 506, 431 

Mordasini, C, Alibert, Y, & Benz, W. 2009, A&A, 501, 1 139 

Morel, P. & Lebreton, Y. 2008, Ap&SS, 316, 61 

Pepe, E, Mayor, M., Galland, E, et al. 2002, A&A, 388, 632 

Perruchot, S., Kohler, D., Bouchy, E, et al. 2008, in Society of Photo-Optical 

Instrumentation Engineers (SPIE) Conference Series, Vol. 7014, Society of 

Photo-Optical Instrumentation Engineers (SPIE) Conference Series 
Pinheiro da Silva, L., Rolland, G, Lapeyrere, V, & Auvergne, M. 2008, 

MNRAS, 384, 1337 
Press, W. H., Teukolsky, S. A., Vetterling, W. T, & Flannery, B. P. 1997, 

Numerical Recipes in C, The Art of Scientific Computing, 2nd edn. 

(Cambridge University Press) 
Queloz, D., Bouchy, E, Moutou, C, et al. 2009, A&A, 506, 303 
Sartoretti, P. & Schneider, J. 1999, A&AS, 134, 553 
Sato, B., Fischer, D. A., Henry, G. W., et al. 2005, ApJ, 633, 465 
Schneider, J. 2003, in SF2A-2003: Semaine de l'Astrophysique Francaise, ed. 

F. Combes, D. Barret, T. Contini, & L. Pagani, 149-+ 
Siess, L. 2006, A&A, 448, 717 
Sing, D. K. 2010, A&A, 510, A260000+ 

Surace, C, Alonso, R., Barge, P., et al. 2008, in Society of Photo-Optical 
Instrumentation Engineers (SPIE) Conference Series, Vol. 7019, Society of 
Photo-Optical Instrumentation Engineers (SPIE) Conference Series 

Valenti, J. A. & Fischer, D. A. 2005, ApJS, 159, 141 

Valenti, J. A. & Piskunov, N. 1996, A&AS, 118, 595 



9 



P. Borde et al.: Transiting exoplanets from the CoRoT space mission 



Institut d'astrophysique spatiale, Universite Paris-Sud 11 & CNRS 
(UMR 8617), Bat. 121, 91405 Orsay, France 
e-mail: pascal . borde@ias . u-psud . f r 

2 Institut d'astrophysique de Paris, Universite Paris 6 & CNRS (UMR 
7095), 98 bd Arago, 75014 Paris, France 

3 Observatoire de Haute-Provence, CNRS & OAMP, 04870 St-Michel 
l'Observatoire, France 

4 Laboratoire d'astrophysique de Marseille, Universite de Provence & 
CNRS (UMR 6110), 38 rue F. Joliot-Curie, 13388 Marseille, France 

5 Institute of Planetary Research, German Aerospace Center, 
Rutherfordstrasse 2, 12489 Berlin, Germany 

6 Laboratoire de l'univers et de ses theories, Observatoire de Paris & 
CNRS (UMR 8102), 5 place Jules Janssen, 92195 Meudon, France 

7 Observatoire de Geneve, Universite de Geneve, 51 Ch. des 
Maillettes, 1290 Sauverny, Switzerland 

8 Department of Physics, Denys Wilkinson Building Keble Road, 
Oxford, OX1 3RH 

9 Instituto de Astroffsica de Canarias, E-38205 La Laguna, Tenerife, 
Spain 

10 Departamento de Astroffsica, Universidad de La Laguna, E-38200 
La Laguna, Tenerife, Spain 

11 Laboratoire d'etudes spatiales et d'instrumentation en astro- 
physique, Observatoire de Paris & CNRS (UMR 8109), 5 place Jules 
Janssen, 92195 Meudon, France 

12 Physikalisches Institut Universitat Bern, Sidlerstrasse 5, 3012 Bern, 
Switzerland 

13 Rheinisches Institut fur Umweltforschung an der Universitat zu 
Koln, Aachener Str. 209, 50931 Koln, Germany 

14 Research and Scientific Support Department, European Space 
Agency, Keplerlaan, NL-2200AG, Noordwijk, The Netherlands 

15 Institute for Astronomy, University of Vienna, Tiirkenschanzstrasse 
17, 1180, Vienna, Austria 

16 Institute of Astronomy, Geophysics and Atmospheric Sciences, 
University of Sao Paulo, Brasil 

17 Thiiringer Landessternwarte, 07778 Tautenburg, Germany 

18 IAG, Universite de Liege, Allee du 6 aout 17, Liege 1, Belgium 

19 Universite de Nice-Sophia Antipolis, CNRS UMR 6202, 
Observatoire de la Cote d'Azur, BP 4229, 06304 Nice Cedex 
4, France 

20 Space Research Institute, Austrian Academy of Sciences, 
Schmiedlstr. 6, Graz, Austria 

21 Wise Observatory, Tel Aviv University, Tel Aviv 69978, Israel 

22 Center for Astronomy and Astrophysics, TU Berlin, 
Hardenbergstrasse 36, 10623 Berlin, Germany 



10 



